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Abstract 



We consider the joint distributions of particle positions for the con- 

■ tinuous time totally asymmetric simple exclusion process (TASEP). 
! They are expressed as Fredholm determinants with a kernel defining a 

o 

00 

■ initial conditions and determine the kernel in the scaling limit. This 
^sO , result has been announced first in a letter by one of us [34] and here 

we provide a self-contained derivation. Connections to last passage 
directed percolation and random matrices are also briefly discussed. 



o 

1 Introduction 



Continuous time TASEP. The totally asymmetric simple exclusion process 
(TASEP) is one of the simplest interacting stochastic particle systems. Its 
particles are on the lattice of integers, with at most one particle at each 
site (exclusion principle). The dynamics of the TASEP is denned as follows. 
Particles jump on the neighboring right site with rate 1 provided that the 
site is empty. This means that jumps are independent of each other and take 
place after an exponential waiting time with mean 1, which is counted from 
the time instant when the right neighbor site is empty. 

On the macroscopic level the particle density, u(x,t), evolves determin- 
istically according to the Burgers equation d t u + d x (u(l — u)) = [32]. 
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Therefore a natural question is to focus on fluctuation properties and on 
large deviations, which exhibit rather unexpected features. The fluctuations 
of the integrated particle current and in the positions of particles are two 
faces of the same coin as to be discussed later. 

The fluctuations of particle positions are sensitive to the initial conditions. 
For example, one can consider particles initially positioned every second site, 
i.e., on TIi. Another possibility would be to consider the stationary measure 
of the same density as initial condition, which is Bernoulli with density 1/2. 
The scaling exponents for particle positions fluctuations are the same for 
the two initial conditions. However, the limiting distribution differ: as we 
will see, in the first case it is the GOE Tracy- Widom of random matrices [41] 
which differs from the stationary case [12]. Thus it is of interest to understand 
which class of initial conditions leads to the same limit distribution. 

The first result in this direction has been obtained for the step initial 
condition. To be precise, denote by Xk{t) the position of particle k at time 
t, where the /c's are integers labelling the particles from right to left. The 
step initial condition is then £fc(0) = —k, k G N. It has been studied by 
Johansson [15] by means of a growth model. In terms of the TASEP, the 
quantity analyzed is the large time asymptotic fluctuations of the position 
of a given particle. For any fixed a G (0,1), the fluctuation of x\ at -\(t) are 
asymptotically governed by the GUE Tracy-Widom distribution, namely 
there are some v = v(a) and b = b(a) such that 

lim P{x [at] {t) < v{a)t + sb{a)t 1/3 ) = F 2 {s). (1.1) 

t — >oo 

The distribution F 2 first appeared in the context of the Gaussian Unitary 
Ensemble (GUE) of random matrices as the distribution of the largest eigen- 
value in the limit of large matrix dimension [40] . 

A natural question is to ask how the positions of different particles are 
correlated, i.e., one considers for fixed but large time t the process k 1— > Xk(t). 
To illustrate known results we focus at k ~ t/4, but the same holds with 
different numerical coefficients for k ~ at, a G (0, 1) 0. (11.11) tells us that 
the fluctuations live on a t 1 ^ 3 scale and it turns out that the position of 
two particles are, on the t 1//3 scale, non-trivially correlated over a distance 
of order t 2//3 . The exponents 1/3 and 2/3 are indeed the ones of the KPZ 
universality class [18], to which the TASEP belongs. Indeed, Johansson [16] 
proves a functional limit theorem in a discrete-time setting. Its continuous- 

1 We choose a = 1/4 because then the term linear in t disappears in (|1.2[) . (|1.3|) . and 
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time version writes 



& = Mu) (1-2) 

where A 2 is known as the Airy process, first discovered in the PNG model 
under droplet growth [29]. 

Besides the step initial condition explained above, at least two other 
situations are of particular interest. One is the stationary initial condition, 
for which the two-point function of the TASEP is analyzed in [12]. The 
second one has deterministic initial conditions leading to a macroscopically 
uniform density, thus called flat initial conditions. The simplest realization 
is to set Xfc(O) = — 2k, k G Z. 

In [34] an important new result has been obtained, making possible the 
asymptotic analysis for such initial conditions. First of all, as expected by 
universality, the position of a particle has fluctuations governed by the GOE 
Tracy- Widom distribution, i*\ [41]. This result is a combination of [11,34] 
and states 

lim P(x m (t) < -st 1 / 3 ) = F 1 (2s). (1.3) 

t — >oo 

More importantly, in [34] the analogue of the Airy process A 2 for flat initial 
conditions, which we denote by A\, is identified. It is a marginal of the 
signed determinantal point process with the extended kernel (11.41) . Here 
signed refers to the non-positiveness of the measure (it does not define a 
probability measure). Explicitly, let B (x,y) = Ai(x + y) and let A be the 
one-dimensional Laplacian; then 

K Fl ( Ul , Sl ;u 2 ,s 2 ) = -(e^-^ A )( Sl ,s 2 )l(u 2 > u r ) + ( e -« lA B oe « aA )( 5l , s 2 ), 

(1.4) 

or, equivalently as shown in Appendix |A] 

K Fl (u u si;u 2 ,s 2 ) = 7= = :exp (-^r — ^2 > «i) 

+Ai(si + s 2 + (u 2 - Mi) 2 ) exp ^(m 2 - u 1 )(s 1 + s 2 ) + ^(u 2 - . (1.5) 

This is in complete analogy with the Airy process A 2 , which is a marginal of 
the determinantal point process defined by the extended Airy kernel. 

In Theorem 12.11 we provide a derivation of the fact that the joint distri- 
bution of particle positions is given by the Fredholm determinant of a kernel. 
This is a general result which is then applied to the flat initial conditions 
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^fe(O) = —2k, k G Z, see Theorem 12. 21 The proper rescaling of particle 
positions is 

xr c (u) = -t- 1/3 (* [t/4+ut2/3] (t) + 2ut^). (1.6) 

With this rescaling, in the limit of large time t, the kernel converges to 
Kf ± as shown in Theorem 12.31 There we show pointwise convergence. In a 
forthcoming paper [3] we will analyze a discrete time version of the TASEP 
and strengthen our result to the convergence of the Fredholm determinants. 
Such a stronger convergence would imply the convergence of (jl.fip to A\ in 
the sense of finite-dimensional distributions. 

As a remark we want to point out that, while periodic initial condition 
does not seem to be accessible by previously known techniques, with the 
new construction both step and periodic initial conditions can be analyzed. 
The technique used so far is the reduction of the model to a determinantal 
point process via the Robinson-Schensted-Knuth correspondence. For further 
references and details on determinantal point processes we refer the reader 
to surveys [14,17,22,37] and the lecture notes [38]. 

Reformulation of the result. The TASEP integrated current at position x 
and time t, J(x, t), is the number of particles which jumped from site x to site 
x + 1 during the time interval [0, t]. Let us label by 1 the right-most particle 
starting at position £i(0) < x. Then P(J(x, t) > s) = P(x s (t) > x + 1). 
Thus the result of this paper translates directly to the integrated currents. 

The TASEP can be mapped to last passage percolation on Z 2 with i.i.d. 
exponentially distributed random variables uj(i,j), i,j G Z. ui(i,j) is the 
waiting time of particle number j to jump from position i — j to i — j + 1. 
There is a slight switch in the point of view. For the TASEP one considers the 
particle positions at fixed time t, while in last passage percolation, one studies 
the last passage time from the origin to points of a given lattice domain 
{i, j G Z 2 ,i + j = t}. These two points of view are closely connected. They 
can be regarded as taking different cross sections in the Bernoulli cone [27]. 
The problem considered in [15] is the point-to-point last passage percolation, 
which corresponds to the step initial condition for the TASEP. Flat initial 
conditions correspond to the point-to-line percolation. 

Finally, the same last passage percolation model can be seen as a one- 
dimensional growth model [15,33], called discrete polynuclear growth model, 
which serves as a discretized model for KPZ growth [18]. KPZ growth is 
discussed in the books [2,23] but for a recent exposition of KPZ universality 
see [27]. 

Universality issue. The TASEP also has discrete time versions. One of 
these is the parallel update rule (see the review [36]) and it is given as follows. 
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At each time step particles jump to the neighboring right site with probability 
p G (0, 1), provided the target site is empty. The jumps occurs independently 
and simultaneously. There are two interesting limits of the discrete-time 
TASEP, namely p — > and p — > 1. The continuous-time TASEP is ob- 
tained by setting the time-unit to p and then take p — > 0. The limit p — > 1 
yields the so-called polynuclear growth (PNG) model, see e.g. Section 2.1.5 
of [9]. There, one has a height function h on a one-dimensional substrate, 
and flat initial condition for the TASEP translates to growth starting with 
h(x,t = 0) = 0, also called fiat PNG. 

By universality the process A\ is expected to appear in the discrete- 
time TASEP and the PNG model as well. Universality has been confirmed 
for step initial conditions and the corresponding PNG model with droplet 
growth [16,29]. Moreover, for flat initial conditions, the limit process should 
be independent on the initial particle spacing. Results in a discrete-time 
version of the TASEP with different initial spacing will be presented in [3]. 

For the flat PNG model, it is known [1,28] that the height at one fixed 
point is GOE Tracy- Widom distributed in the limit of large time t. Thus, 
on the basis of the result for the TASEP with alternating initial conditions, 
see Theorem 12.31 we can conjecture the behavior of the flat PNG model. 

Conjecture 1. The properly rescaled height function of the PNG model with 
flat initial conditions converges, in the large time limit, to the process A\ . 

The scaling exponents are the same and the coefficients can be determined 
by matching with the PNG droplet case. 

Finally, let us discuss the connection to random matrices. For the TASEP 
with the step initial condition, the one-point asymptotic distribution is the 
GUE Tracy- Widom distribution, F 2 , and the whole limit process is the Airy 
process A 2 - The derivation uses an extension of the model to a multi-layer 
version. The Airy process arises also for a GUE matrix diffusion, the so-called 
(3 = 2 Dyson's Brownian Motion [6]. The motion of the properly rescaled 
largest eigenvalue converges to the Airy process. The connection extends to 
finitely many of the largest eigenvalues which have the same limiting behavior 
as the first top layers in the multi-layer PNG model. 

For the TASEP with flat initial condition, the one-point distribution is 
the GOE Tracy- Widom distribution and the limit process is A\. At this 
point it is tempting to conjecture that the evolution of the largest eigenvalue 
of a matrix which follows (3 = 1 Dyson's Brownian Motion has the same 
limiting behavior as the surface height for flat PNG, namely the A\ process. 
The correspondence at the level of top eigenvalues for GOE and the top 
layers of the multi-layer flat PNG at a fixed position has been proven in [8] , 
making the conjecture even more plausible. Knowing the analogue of the 
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Airy process for random growth with flat initial conditions one can guess the 
result for (3 = 1 Dyson's Brownian Motion [6]. 

Conjecture 2. The evolution of the largest eigenvalue of N x N matrices 
for (3 = 1 Dyson's Brownian Motion converges, in the limit N — ► oo and 
properly rescaled, to the process A\ . 

Again, the prefactors for the scaling can be easily calculated by matching 
the known one-point distributions and the behavior of joint distributions at 
short distances. This conjecture concerns only the largest eigenvalue and with 
this approach we are unable to make a conjecture for the other eigenvalues. 

To make Conjecture [2] more transparent, we explain it in the simpler case 
of the two-matrix model. There, one considers two real symmetric N x N 
matrices, M(0) and M(t), with joint distribution 

_ ^_j^ y M{0)AM{t) (17) 

where q = exp(-t/2N) and dM(-) = Yli<i<j<N dM(-) id . Let A max (0) and 
A max (t) be the largest eigenvalues of M(0) and M(t). These eigenvalues 
fluctuate on a scale of order iV 1//3 and are non-trivially correlated if one 
chooses t ~ iV 2 / 3 . Then Conjecture [2] means that, properly rescaled, the 
joint distribution of A max (0) and A max (t) converges to the two-point joint 
distribution of the process A\ in the N —>■ oo limit. 

We also refer to the surveys on the question of universality in mathemat- 
ics and physics [5] and on connections between different models, including 
random matrices [10]. 

Outline. The paper is organized as follows: In Section [5] we state the 
main result. In Section [3] the kernel of the signed determinantal point pro- 
cess describing the joint particle distributions is derived. The kernel involves 
an orthogonalization which is carried out in Section H] for the case of al- 
ternating initial conditions. In Section [5] we prove the convergence of the 
properly rescaled kernel to the kernel Kp 1 . In Appendix [X] we explain how 
the compact form of the kernel is derived, and in Appendix [B] we explain how 
the orthogonalization can be carried out using classical Charlier orthogonal 
polynomials. 
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2 Model and results 

In this paper we consider the continuous-time totally asymmetric simple 
exclusion process (TASEP) on Z. At any given time t, every site j G Z can 
be occupied at most by one particle. Thus a configuration of the TASEP can 
be described by rj = {rjj,j G Zi\r)j G {0, 1}} G Q. rjj is called the occupation 
variable of site j, which is defined by rjj = 1 if site j is occupied and rjj = 
if site j is empty. 

The dynamics of the TASEP is defined as follows. Particles jumps on 
the neighboring right site with rate 1 provided that the site is empty. This 
means that jumps are independent of each other and are performed after an 
exponential waiting time with mean 1, which starts from the time instant 
when the right neighbor site is empty. More precisely, let /: Q — > R be 
a function depending only on a finite number of rj/s. Then the backward 
generator of the TASEP is given by 

Lf(v) = 5>( X " V j+ i)(f(rf> j+1 ) - MY (2-1) 

Here rf^ +1 denotes the configuration r\ with the occupations at sites j and 
j + 1 interchanged. The semigroup e Lt is well-defined as acting on bounded 
and continuous functions on Q. e Lt is the transition probability of the 
TASEP [21]. 

Joint distributions 

Let us start at time t — with N particles at positions Un < ■ ■ ■ < V2 < Hi- 
Then the main result is the joint distribution of any subset of these particles 
at time t > 0. It turns out that it can be described by a signed determinantal 
point process, where signed refers to the non-positiveness of the measure. 

Theorem 2.1. Let er(l) < er(2) < . . . < a(m) be the indices of m out of the 
N particles. The joint distribution of their positions x a ^(t) is given by 

m 

P ( fl {^WW - = det ( 1 ~ XaKtXa)t?({<TQ.),...,*(m)}xT) ( 2 -2) 

k=l 

where Xa(o~(k),x) = t(x < a k ). K t is the extended kernel with entries 

K t {n u x^x 2 ) = -^\ Xl ,x 2 ) + Kl-n 2+l ^iM 2 (^) (2.3) 

i=0 
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where 

*"■"-'<»„*)= (:;:::!). <^> 

*?<*> = ^ / 4^^^-". (2-5) 

27T1 J ro W l+i W x 

and the functions <&™(x), i = 0, . . . , n — 1, form a family of polynomials of 
degree < n satisfying 

£ *?(*)<&?(*) = «5 W . (2.6) 

sea 

T/ie pai/i T in the definition of is any simple loop, anticlockwise oriented, 
which includes the pole at w = but not the one at w = 1 . 

The dependence on the set {y{\ is hidden in the definition of the $™'s 
and the \&"'s but is omitted, since the set is fixed in the following. 



Alternating initial configuration 

Now we consider alternating initial configuration, namely 

f 1, if i is even, _ N 

»<(°)={o; it sis odd. < 27 > 

The alternating initial configuration can be obtained by taking 2N particles 
around the origin, for example at positions 27L H [— 2N, 2N — 2], and then 
taking the limit N — > oo. In Lemma [4. II we do the orthogonalization, i.e., 
construct $"'s which satisfy (12.61) for this special case, from which the kernel 
K t is obtained. 

Theorem 2.2. Let particle with label rij start at —2i, i G Z. At time t, 
the particles are at positions X{. The kernel $2.31) for the alternating initial 
configuration is given by 

,'x 1 -x 2 -l\ -1 / , (1 + v) xa+ni+na tn , 2v] 
K t (n u x x - n 2 , x 2 ) = - J + ^ * dv ) '^ e-« 1+2v > 



n 2 -n 1 -\J 2ni J To [- v )^+ni+n 2 +i 

^ _ "" (2-8) 

where T is any simple loop, anticlockwise oriented, which includes the pole 
at v = but do not include v = — 1. 
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Scaling limit 

The particle density is 1/2 and since particles jump to the right with rate 
1 provided the site is empty, the mean speed of the particles is 1/2. Let us 
number the particles from right to left with yi(0) = as reference point, 
i.e., yi(0) = —2(i — 1), % G Z. Then the particles which at time t are close 
to x = are the particles with numbers close to t/4. From universality we 
know also that the scaling exponent for fluctuations should be 1/3 and the 
one for spatial correlations should be 2/3. Therefore, the scaling limit to be 
considered is 

Xi = -2u t t 2/3 - Sit 1/3 , m = t/4 + mt 2/3 . (2.9) 

Remark: The scaling exponents for this model are determined by the re- 
quirement that there is a non-trivial limit. The numerical factor in front of 
t 1 / 3 is chosen so that the single-time kernel has a simple form, Ai(x + y). 
The numerical factors for the t 2//3 terms are set in such a way that the prop- 
agator in (11.41) is generated by the Laplacian without additional prefactors. 
Universality argument is not needed to obtain the result, but it is useful to 
predict the correct answer. 

In Section [3] we carry out the asymptotic analysis for the pointwise con- 
vergence of the kernel, with the following result. 

Theorem 2.3 (Pointwise convergence of the kernel). Let Xi,ni,x 2 ,n 2 be 
reseated as in K2.9\) . Then, for any s\, s 2 , U\, u 2 G R fixed, 

\im K t (n 1 ,x 1 ;n 2 ,x 2 )t 1 ^2 X2 -^ = K Fl ( Ul ,s i; u 2 ,s 2 ) (2.10) 

t— »oo 

where the extended kernel Kp 1 is given in ( ti.5]) . 

In this paper we do not perform the asymptotic analysis necessary to get 
convergence of the Fredholm determinants. We will do the complete analysis 
in a discrete-time version of the TASEP in a forthcoming paper [3], from 
which the continuous time limit follows as a corollary. Nevertheless, it is 
instructive for the reader to see the implications of the convergence of the 
Fredholm determinant. 

Let Ai be the process with m-point joint distributions at 
u\ < u 2 < . . . < u m given by 

m 

P( f)Ml(«fc) ^ = det ( 1 ~ X-^F 1 X-)za({u 1 ,...,« ra }xR) (2-11) 

fc=l 
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where Xs{uk,x) = l(x > s^). The convergence of Fredholm determinant 
would then imply 

hm = Mu), (2.12) 

with the convergence understood in the sense of finite-dimensional distribu- 
tions. 



3 Signed determinant al point process 

In this section we prove Theorem 12.11 Consider the TASEP with iV particles 
starting at time t = at positions < . . . < y 2 < y\- The first step 
is to obtain the probability that at time t these particles are at positions 
xjsr < . . . < X2 < Xi, which we denote by 

G(xi, . . . ,x N ;t) = P((x N ,...,x 1 ;t)\(yN,---,yi;0)). (3.1) 

This function has been determined before using Bethe-Ansatz method [35]. 

Lemma 3.1 (Schutz [35]). The transition probability has a determinantal 
form 

G(xt, ...,x N ;t) = det(Fj_ i (a; A r + i_j - y N +i-j, t))i<i,j<N (3.2) 



with 



Fn{x , t) = t^I ^^e'<->, (3.3) 



27ri J Tq 1 w w x ~ v 



where r 0> i is any simple loop oriented anticlockwise which includes w = 
and w = 1. 

This representation of the transition probability was utilized to study the 
current fluctuations in [24,30]. To study the joint distribution, we need a de- 
composition of G(x\, . . . , xn', t) given in the next lemma. This decomposition 
is obtained using only the recurrence relation 

F n - X (x, t) = F n (x, t) - F n (x + 1, t) (3.4) 

and its integrated form 

F n+1 (x,t)=J2 F n(y,t). (3.5) 

y>x 
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Figure 1: Graphical representation of the domain of integration T> for iV = 4. 
One has to "integrate" out the variables x\, i > 2 (i.e., the black dots). The 
positions of x\, k = 1, . . . , N are given (i.e., the white dots). 



Actually, (13.51) comes from (13.41) and the fact that lim^^ F n (y, t) = fast 
enough. The other property needed to obtain Theorem 12.11 is the following. 
Comparing (12.51) and (13. 3p . we have 

^(x) = (-l) k F_ k (x-y N ^ k ,t) (3.6) 

for k > 0. Notice that, for n = —k < 0, (13. 3p defining F n has actually only 
one pole at w — 0. We then get the relation 

F n+1 (x,t) = -J2 F n(y,t), (3.7) 

y<x 

which translates into 

<- k (x) = J2 xS, N + +l^- ( 3 - 8 ) 

y<x 

In the definition of the ^j^'s in Theorem 12 .1} the path T includes only the 
pole at the origin, exactly because we need (13. 8j) to hold also for k < 0. 

Lemma 3.2. Let us denote x k = x\, k = 1, . . . , N. Then 

G(xi, . . . ,x N ;t) = y j det(F_ i (x^ 1 - y N -j, Q)o<i,j<jy-i (3.9) 
v 

where the sum is over the following set 

V = {xj, 2<i<j< N\xi > xi+\ x\ > xill}. (3.10) 
See Figured for a graphical representation ofV. 
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This lemma is actually more general as shown in [34]. By applying just 
the recursion relation (I3.5P the domain of summation in Lemma I3~2l would be 
V = {x{, 2 < i < j < N\x{ > xjZi} instead of D. The reduction of the sum- 
mation domain to T> uses only the antisymmetry of the determinant. Thus 
the same holds for any antisymmetric function /. It might be interesting in 
other applications, so we state it explicitly. 

Lemma 3.3. Let f an antisymmetric function of{xi, . . . , x^}. Then, when- 
ever f has enough decay to make the sums finite, 

^/(xf,...,x^) = ^/(xf,...,4) (3.11) 



V V 



where 



V = {.'••;. 2 < i< j< N\x{ >xy\rj>xj .'}. 

V = {x{, 2 < % < j < N\xj > x{zl}, (3.12) 

and the positions x\ > x\ > . . . > Xi being fixed. 

Proof of Lemma \3.2l The proof consists in applying the property (13.51) itera- 
tively and using the multilinearity of the determinants. From Lemma T3. II we 
have 



G(x{,...,x?;t) = det 



F (xf-y N ,t) ■■■ F_ N+ i(x? - y x ,t) 

F N -i(x{-y N ,t) ••• F (x\-yi,t) 
The first step is to rewrite the last row as 

[ F N . 2 (x 2 2 -y N ,t) ■■• F^(x 2 2 -y u t) ]. (3.14) 



(3.13) 



The second step is to apply the same procedure to the last and second to 
last rows, which become 

Yl E I F N-3(4-VN,t) ■■■ F_ 2 (x 3 3 -y u t)] (3.15) 

and 

[F N -3(x 3 2 -y N ,t) ■■■ F„ 2 {x 3 2 -yi,t)] (3.16) 
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respectively. At this point we have 

F (x?-y N ,t) 



= E E E det 



F Ar _ 3 (a;| -y N ,t) 
F N _ 3 (xl -y N ,t) 



F_ N+1 (xf -y u t) 



F- 2 {x{ 
F-2{x\ 



yi,t) 
yi,t) 
yi,t) 



(3-17) 

The determinant is antisymmetric in the variables (xf, xjj), therefore the con- 
tribution of the symmetric part of the summation domain of E/x§>s§ ^2 x 3 >x j 
is zero. Since ^ x 3> x 2 ^ x s> x 2 = Ei'^ EiSe[ij,^) + E^>^> tne 

symmetric part of the domain is {xjj > x\,x\ > x^}, thus the contribution 
coming from the last sum is zero. 

We iterate the same procedure. More precisely, for k = 4, . . . , N, we 
apply (13. 5p to the last (k — 1) rows. The new summing variable for the last 
row is denoted by x%, the second last row x^_ v and so on. Finally, we can 
delete the sums over the symmetric domain in (x§, . . . , x|). In this way we 
get the result 



G(x\,...,x?;t) = J2 det 



v 



F (x?-y N ,t) 



F^ N+1 (x? -y h t) 



F (x% 



VN,t) 



F_ N+1 (x% -yi,t) 



(3.18) 
□ 



This is the decomposition used in [34]. The integrations variables 
{xf,i = l,...,n} can be interpreted as the positions of particles labelled 
by i — 1, . . . , n at time n. For example, the trajectory of par- 

ticle 1, see also Figure EJ This is just a mathematical construction which 
should not to be confused with the real TASEP particles and the natural 
time in the TASEP positions, which at this stage is just the fixed parameter 
t. At time n there are n particles at positions x™, . . . , x™. At time n + 1, they 
jump to a randomly uniformly chosen position satisfying x£ +1 G [x^_ x ,x^) 
with the (n + l)st particle added at position x™+}(> x"). Then the weight 
of a configurations of x"'s is given by 

W({x^;l <i<n<N}) 



(3.19) 



N 

( JJdet(l(xf _1 > x"))i<M<n) det(F_ j (x^_ 1 - y N -j,t))o<i,j<N-i, 



n=2 



where we set x" 1 = oo. The products of determinants in (13.191) might look 
complicated. However, one can verify that whenever some of the Xj's do not 
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Va V3 2/2 2/1 



N + l 




N = 4 
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1 
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x 



Figure 2: LGV scheme for N = 4. The LGV graph is left/up- left directed 
with weight 1 for each edge. From iV to N + 1 the transitions are the F's. 

satisfy the inequalities of T>, then at least one of the determinant vanishes. 
On the other hand, if the set of x^s belongs to T>, then each determinant has 
value 1. 

The form of the weight suggests that the correlation functions could 
be determinantal. It is like to have a Lindstrom-Gessel-Viennot (LGV) 
scheme [42] , see [39] for a nice exposition, with a sort of reservoir of particles 
at oo and at each time-step a new particle is introduced. The LGV scheme is 
a sort of generalization on a class of directed graphs of the Karlin-McGregor 
result for diffusions [19]. Determinantal form of correlation functions ap- 
peared in different contexts [4,7,13,16,25]. Although we do not use the LGV 
scheme in the proof, it might be interesting for the reader to see how the 
weight (I3.19P can be described in this framework. The situation given by 
the weight (13.191) corresponds to the limit U — ► oo of the system with fixed 
number of particles illustrated in Figure EJ As U — > oo, the extra particles 
are not seen, they goes to oo in a sort of reservoir. 

The proof of Theorem 12. II is an application of the following Lemma, which 
is proven by using the framework of [4]. 

Lemma 3.4. Assume we have a signed measure on {x™,n = 1, . . ., N, % = 
1, . . . , n} given in the form, 



where x™ +1 are some "virtual" variables and is a normalization constant. 
If Zn 7^ 0, then the correlation functions are determinantal. 

To write down the kernel we need to introduce some notations. Define 



1 



Y[ det#„(. 



X™, < +1 )]i<M< n+ i detftfj^ay )]i<i,j<N, 



(3.20) 




* (f> n2 -i)(x,y), n x < n 2 , 
ni > n 2 , 



(3.21) 
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where (a * b)(x, y) = XLez a ( x ' z )H z i v)> an d, f or 1 < n < N , 

^(x):=(r' N *K- 3 )(y), J = 1,2,..., iV. (3.22) 

Set (/>()(£?, x) = 1. Then the functions 

{(0o * (1,n) )« x), . . . , (0 n _ 2 * (n - 1 ' n) )(x^, x), </> n _i«-\ x)} (3.23) 

are linearly independent and generate the n- dimensional space V n . Define a 
set of functions j — 0, . . . ,n — 1} spanning V n defined by the orthog- 

onality relations 

^n^](x) = 5 hJ (3.24) 

IE 

/or < i,j < n — 1. 

Under Assumption (A): n (x^_|_^, x) = c n $o"' +1) (x), /or some c n 7^ ; 
n = 1, . . . , AT — 1, t/ie kernel takes the simple form 

"2 

K( ni , Xl ; n 2 , x 2 ) = -<f>^\ Xl , x 2 ) + ^ Kl-kW^Z-M- ( 3 - 25 ) 

fc=i 

Remarks: Without Assumption (A), the correlations functions are still 
determinantal but the formula is modified as follows. Let M be the N x N 
dimensional matrix defined by [M]i t j = ((fii-i * <p l,iV ) * 1 ^^_j)(x t i ~ 1 ). Then 

K(n 1 ,xi;n 2 ,x 2 ) (3.26) 

n 2 N 

-^\x h x 2 ) + J2Kl- k (xi) YslM^Mh-i * ^ l ' n2) M-\x 2 ). 

k=l 1=1 

The analogue of the determinantal representation (13.21) for particle- 
dependent hopping rates has been recently obtained [31]. Lemma [3.41 might 
be applied in this context too. 

Proof of Lemma We apply Proposition 1.2 of [4] and we try to stick 
as much as possible to the notations therein. Let, for n = 1,...,N, X^ 
denote the space of {x?, i = 1, . . . , n}, 2) = X^ U . . . U X^ N \ and let 
X = {xi, x\, . . . , xj^ -1 } U 2) be the space on which our measure (13.201) is de- 
fined. Let X (ra ' n+1 ) be the matrix with entries 

[T^] hj = 0("' n+1 )(x^,x™ +1 ), l<i,j<n+l (3.27) 

and 

[* W )<j = *N-j{*?)> l<iJ<N. (3.28) 
15 



Then the weight (13.201) is proportional to the determinant of 

-T {1 ^ ■•■ 
-T( 2 > 3 ) ••• 

















_rp(N-l,N) 





(3.29) 



We are interested in the measure on 2J only, thus we change the ordering 
by putting the variables x\, . . . , x^" 1 at the beginning. Let us define the 
n x (n + 1) matrix l^L n+1 ) by 

[W [n , n+1) ] i7j = ( "' n+1) (xr, x] +1 ) } 1 < i < n, 1 < j < n + 1, (3.30) 
and the N x (m + 1) matrix E m by 



[E, 



m\i,j 



x™^,x? +1 ), i = m + 1,1 < j < m + 1, 



0, otherwise. 



(3.31) 



Then the weight (13.201) is proportional to a suitable symmetric minor of L, 
with 





' 


Eo 


E 1 


E 2 


En-i 










-W[l,2) 








L = 






























—W[N-1,N) 



















(3.32) 



By Proposition 1.2 of [4], the point-measure on 2) is determinantal with 
correlation kernel given by 



2)x2) 



(3.33) 



provided that the partition function Zn ^ 0. With the decomposition of 
X = {xi, x\, . . . , £jv _1 } U 2), we have a block decomposition of L as 



B 
C D 



(3.34) 



with B = [E , . . . , E N -i], C = [0, . . . , 0, ^W]*, and D equal to L without 
the first line and column of the block representation (13.321) . Let D = 1 + Do, 
then (see, e.g., Lemma 1.5 of [4]) the kernel is given by 



K = t — D 



-i 



D^CM-'BD 



-1 c n-1 



M = BD- X C 



(3.35) 
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D 1 was already computed in Lemma 1.5 of [4], with the result 

1 W [1]2) ■■■ W [ltN) 

oi'-. ; 



D 



-i 



1 



(3.36) 



where 



W [n , n+ i) ■ ■ ■ W[ m - 1>m ), m> n, 
"•'"'-1 0, m<n. 

Thus the (n, m) block of 1 — D^ 1 is W\ n>m ). Next, we have 



(3.37) 



D~ C 



(3.38) 



and 



BD- 1 = [ E E Q W m + Ei ■■■ Y,k=i E k -iW [KN) + E N ^ } . (3.39) 
Therefore the (n, m) block of the correlation kernel is given by 



m—l 



J m—1 



(3.40) 



K {n > m) = -W [n . m) + W^n^M- 1 ( E k-iW[k,m) + E n 

fc=i 

ag (E2U) one gets [W [ritm) ] id = 0( n > m )(x™, xf). Moreover by (E22D we 
e 



It remains to evaluate the last part of (13.401) . For the following N xm matrix 
we have 



m— 1 



fc=l 







ii-i * (i ' m) )(x*- 1 ,x J m ), 1 < i < m, 

m + 1 < z < N. 



(3.42) 

Notice that the functions in (13.421) form a basis of V m . Thus we can define 
a m x m matrix B m which does a change of basis to {$™_ x (x), . . . , $™(x)}, 
namely 



W ftm) )(4V) = £[iU^ 



(3.43) 



z=i 
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We multiply this equation by ^2 x ^m-j( x ) anc ^ OD t arn 

[s m ] iJ = (^_ 1 *^*^_ j )(zr 1 ). 

In particular, we have = M. Let us define the N x m matrix 



^{xf), l<i<m, 
0, 



m + 1 < % < N. 



(3.44) 



(3.45) 



Then 



m—l 



fc=i 



5 m 




(3.46) 



Assume the condition (B): (13.461) = for m = 1, . . . , N. Then we get 
the simple form of the kernel, (13.251) . of the Lemma. However, this is not 
always the case. For 1 < i,j < rn, we obtain using (I3.22p . 



[B m ]ij = (0,-i * (i ' m) * VZ-jKJr 1 ) = (0,-i * <P {l ' N) * K-j)^ 1 ) = [B N k r 

(3.47) 

for some (N — m) x m matrix Q m . 



Thus we can write B 



N 



B r , 



Qm * 

By multiplying on both sides by Bn the condition (B), we see that (B) is 
equivalent to 



k=l 



<f> m 

m\i,k^ m -k 







(3.48) 



for all x and for all % = 1, . . . , N — m, and for all m = 1, . . . , N. But the 
functions $™_ fc (o;) f° rm a basis of V m , thus (B) is fulfilled iff Q m = for 
all m = 1, . . . , N. Thus (B) is equivalent to the condition B^ is an upper- 
diagonal matrix. 

Assume B m upper diagonal for some m. This is verified for m = 1 where 
Bi = I. Then by 03.471) [B rn+1 ] Lj = [B m } itj for % = 1, . . . , m, and B m+l is still 
upper-diagonal iff 



[B 



m+ljm+l ,j 



(3.49) 



c m 7^ because ^ 0. Finally, the orthogonal relations (13.241) imply that 
m (x™ +1 ,x) = c m $™ +1 (x), which is Assumption (A) of the Lemma. □ 

Proof of Theorem \2.1[ It consists in an application of Lemma 13.41 with 

4> n {a%,x? 1 ) = l($>xp 1 ), n = l,...,N-l, (3.50) 
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and 

K-ifr) = (-l) N ~ l F_ N+l (x - Vi ,t), i = 1, . . . , N. (3.51) 

An important point is that in (13.31) the functions F^s are defined by an 
integral enclosing w = and w — 1. At this stage, we have only functions 
FSs for i < 0. In this mentioned around (13.6j) . w = 1 is actually 

not a pole, thus the weight (I3.19P and the weight (I3.20p with the above 
replacements are proportional. The definition of the ^'s using only one pole 
fit exactly in the framework of Lemma 13.41 In fact, by (13.81) . we have the 
composition rule 

(0**£i-iX x ) = *n-i( x )> (3.52) 

which gives (13.221) by iterations. In our setting, if we sum up all the variables 
{x™, 1 < m < n,l < j < m}, we get a Vandermonde determinant in the vari- 
ables x™'s. Thus the space V n of Lemma I3~4l is generated by {1, x, . . . , x n_1 } 
and are polynomials of order at most n — 1. A simple computation using 
(1231) leads to 

W"(x) = {°' ^ = 1 >--->™- 1 > (3.53) 

which, together with (I3.24p leads to $q(x) = 1 = 0„_i(oo,x). Thus we have 
a determinantal system with kernel (I3.25j) . which can be rewritten as 

n 2 — 1 

K t (n 1 ,x 1 ;n 2 ,x 2 ) = -^\x u x 2 ) + ^ Kl- m+ i(^M 2 (^) ■ (3.54) 

i=0 

Since in this paper we explain the detail of the derivation in [34] , it is use- 
ful to point out a difference. There one does not obtain directly 0( ni>n2 ) as in 
Theorem 12.11 There the one-time transition becomes <p n (x,y) = —t(y > x) 
and the representation (I3.3P has to be used instead. The final form of the 
kernel (eq. (13) in [34]) comes from splitting the contribution from the pole 
at w — 1 and the remainder. In the geometric picture, it corresponds to 
have the conjugate LGV graph with reservoir of particles at — oo instead of 
at +oo. □ 



4 Orthogonalization 

In order to get the kernel for the alternating initial configuration, i.e., for 
the case where particles initially occupy exactly the sublattice 2Z, we start 
with a finite number of particles, 2N. In the second step we will focus on the 
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Figure 3: Trajectories of the 2N particles. The black dots are the initial 
positions and the white dots are the positions of the particles at some later 
time t. This is a scheme leading, in the N — > oo limit, the alternating initial 
configuration on Z. 



region where the iV-dependence vanishes. In this way we will get the kernel 
for the system we actually are interested in. 

Consider the case where at time t — there are 2N particles placed every 
second site centered around the origin, see Figure El namely 

Vi = 2N-2i, t = l,...,2N. (4.1) 

From Theorem 12.11 the kernel is known once the orthogonalization is carried 
out. Here we state the result and a short proof. In Appendix [B] we explain 
a constructive way of obtaining needed functions using Charlier orthogonal 
polynomials. 

Lemma 4.1. The functions ^ k (x) an< ^ ^k( x ) have the following integral 
representations. Let z = x + 2n — 2N. Then 

W( x ) = 1 L A ™ e (v-i)tf{ w - \) w f (4.2) 

feV ' 27T1 J Vo W z+1 U ' ' V ' 

and 

27ri % v (v(i + v)f [A - 6) 

where Tq is an anticlockwise simple loop enclosing only the pole at 0. 
Proof. We have 

= (-l) k F_ k (x - y n _ k , t) = (-l) k F_ k (z - 2k, t). (4.4) 
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Then (13.31) leads directly to (14. 2p . Next we prove that (14. 3 p satisfy the or- 
thogonality relation (I3.24p . Since ^(x) = for z(x) < 0, we have 

y^ n k (x(z))^(x(z)) = / dw^-^dw - l)w) k (4.5) 

T^n 2711 Jr n 



z>0 



{-iy r dv (i + 2t») ^ (v + 1) 



27ri J Tq v e vt (v(v + 1)) 



— Y 



z-1 



W 

z>0 



z+l 



provided that the integration domain satisfy |1 + v\ < \w\. The last sum 
gives 

^(v + iy- 1 i . . 

> = . (4.6) 

w z+1 (w - (1 + v))(l + v) v ; 

Thus ( 14 .5p has a simple pole at w — 1 + v, and once the integral over w is 
computed, we get 

f_i \ k +j r 1-1-97) 

E ] (*(*)) = V— f ^itrz^^ 1 + ^ 4 - 7 ) 

z>0 1/10 v y 

The final step is a change of variable. Let u = v(l + v). Then 

du = (1 + 2w)dt» (4.8) 
and the integral is again around 0. Thus we get 

E nwmw*)) = t ^ f = ( 4 -9) 



2>0 



□ 



Once the orthogonalization is made, we determine the kernel of Theo- 
rem ES 

Proof of Theorem \2.2l We need to derive the formula for the first term (the 
main part) of the kernel. For convenience, we first shift the integrating 
variable of ^ to go around —1 by setting u — w — 1. This leads to 

with z = x + 2n — 2N . We start with particles at positions yi = 2N — 2i, 
i = 1, . . . , 2N. The main term in the kernel writes, with Zi = x\ + 2(raj — N), 

— 1 r_i)»i-"a r (I + 2v){l + v) Z2 



e vt (v(l + v)) n2 



X 4> du / V r ^ ; r ; r (4.11) 

Jr^ (l + u)^ +1 u(l + u)-v(l + v) V ; 
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provided that the integration paths satisfies (a) \u(l + u)\ > \v(l + v)\ and 
(b) u — is not inside the contour T_i. To obtain this expression we first 
take the finite sum over k inside the integrals, and secondly we extend it to 
k = — oo. This can be done since the sum is absolutely summable because of 
(a) and we do not create new poles inside the integration contours because 
of (b). For example, we can set Fi by |1 + u\ = 1/2 and take T to be a 
contour with \v\ small enough. 

To obtain the kernel for the alternating initial configuration we focus 
on the Xj's far enough from the right-most particle so that the system in 
the considered region becomes independent of the fact that we have only 
a finite number of particles. This is obtained when Z{ < i.e., whenever 
u = — 1 is not anymore a pole. This condition is satisfied for any fixed Xj 
(i.e., around the origin) and any finite time t by taking iV large enough. In 
fact, it corresponds to taking n« — iV = 0(1) in N. In this case we are left 
with one simple pole at u = — 1 — v. Denote rij = N + uii, then Z{ = Xj + 2rrii 
and the main part of the kernel becomes, for any Xj's as iV — > oo, 

n _ 2 ~} _1 C (1 _|_ „,\x 2 +m 1 +m 2 

Finally, by relabelling the particles we obtain (12.81) . □ 



5 Asymptotic analysis 

In this section we do the asymptotic analysis for the alternating initial con- 
ditions and prove Theorem 12.31 Just to remind, the scaling limit we have to 
consider is 

m = t/4 + Ul t 2/3 . (5.1) 

Proof of of Theorem \2.3i The pointwise limit of the first term is quite easy 
to obtain. Let us set a = (u 2 — ui)t 2 ^ 3 — 1, b — (s2 — siji 1 / 3 + 1, and e — b/a. 
Then we have to compute 

* i/3 ( a(2 ; £) ). (5.2) 

We simply use x\ = \/2nx exp(x ln(a;) — a;)(l + 0(x~ 1 )). Since for s\, s 2 is a 
bounded set, e — >• as t — » oo, we have that 

t i/ 3 ^(2 + e)\ = tl/32Xl _ X2 _l exp( _ 6 2 /4a)(1 + Q{e)) (53) 
\ a J ^J^Txa 
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and by replacing back the expressions of a and b we get, 



lim t 1/3 



t^oo 



X\ — X 2 — 1 

n 2 — ni — 1 



1 



a/47t(m 2 - «r 



: exp 



S2 - Si 



4(w 2 - Mi) 



(5.4) 



Next we analyze the second of the kernel (12. 8p multiplied by t 1//3 . This 
writes 



_ t l/3 

2vri 



r 



d^exp (t/ ( V ) + t 2 / 3 /^) + t^f 2 {v) + / 3 (v)) (5.5) 



with 



fo(v) 
fi(v) 

h(v) 

fs(v) 



2 h i—)- i - 2v - 

-(u 2 -ux)ln((l + v)(-v)) 
—s 2 ln(l + v ) + Si ln(-f ), 
— ln(— 1>). 



(5.6) 



To do a steep descent analysis we first have to find the stationary points of 
fo{v). Simple computations lead to 



d/o(tO 



[l + 2vf 



dv 2v(l + v) 

which has a double zero at v = —1/2. Moreover 
d 2 fo(v) n d 3 / (^) 



(5.7) 



dv 2 



v=-l/2 



0. 



dt> 3 



v=-l/2 



16. 



(5i 



The steep descent path r used for the analysis, shown in Figure HI is 
given by T = V T 2 V T 3 with 



r 1 

J- n 



■i + ice-^ice [0,1/2]}, 



r 2 

1 



{«=- 

{v = -\e ie ,6e [7r/3,57r/3]}, 
^ = {v = -I + (1/2 -w)e^\we {0,1/2}}. 

Let us verify that T is actually a steep descent patrH. On Tq 

dRe(/ o )(0) 4 sin 0(1 



(5.9) 



o- 



cos( 



dO 



5 - 4 cos 



(5.10) 



2 For an integral / = j^dze*^^, we say that 7 is a steep descent path if (1) Re(/(z)) 
is maximum at some z S 7: Re(/(z)) < Re(/(zo)) for z e 7 \ {z } and (2) Re(f(z)) is 
monotone along 7 except, if 7 is closed, at a single point where Re(/) reaches its minimum. 
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Figure 4: The steep descent path T used in the asymptotic analysis. 



Therefore the real part of f is stationary only at 9 = 0, ir and the maximum 
is at 9 = 0, the minimum at 9 = it. By symmetry we need to check only on 
rj. We find 

dRe(/ )H 8 W 2 (1 + 2 W 2 ) 

dw (1 + 2w + Aw 2 ){\ -2w + Aw 2 ) 1 ' ' 

which is strictly negative except at w = where is zero. Thus r is a steep 
descent path. 

Consider now the piece of the path 1^ = {z G r°||z + 1/2 1 < 5}. Let us 
denote F(v) = exp(tf (v) + t 2/3 /i(^) + t 1/3 f2(v) + f 3 (v)). Then since T is a 
steep descent path, 



27ri 



+1/3 r +1/3 r +1/3 r 

J) dvF(v) = / dvF(v) + / dvF(v) 

Jt 2m Jr* 2tti Jr \r* 

t l/3 



[ dvF{v) + F{-l/2)0{e-^) (5.12) 
Jt* 



2ni 

for some pi > (in our case, fx ~ 5 3 ). The precise expression of F(— 1/2) is 
_p(_ 1/2) — 2 2 ( M 2-«i)t 2/3 +(s2-si)t 1/3 +i _ 2 a; i- :E 2+i ^5 ^ 

For the integral on Tq we can apply Taylor development. On Tq we have 
v — — 1/2 + e ±t7r / 3 w, < u> < 5, and we obtain 

/„(«) = -^ 3 + 0(w 4 ), 

/i(u) = 2( M2 - Ml )ln2 + 4( M2 - Ml )e 27ri/3 w; 2 + C(7i; 3 ), 

/ 2 (u) = (s 2 -Si)ln2-2( Sl + s 2 )e i7r / 3 7i; + 0(7i; 2 ), 

/ 3 (u) = In 2 + OH (5.14) 
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where the error terms 0(- ■ ■) are uniform for Sx,S2 in a bounded set. Set 
fi(v) to be the expressions fi(v) without the error terms, and similarly F(v): 
F{v) = exp (tf (v) + el*f x {v) + t^f 2 {v) + f 3 {v)). Then 



dvF(v) = —— dvF(v) + — — I dv(F(v) -F{v)). (5.15) 



27ri l r s 2ix'\ l r s 2tti 

To estimate the second integral, we use the inequality \e x — 1| < e' z '|x|. Thus 



_ t l/3 



t !/3 



d«(F(v) - F(v)) < / dw\F(v(w) = e~' m/3 w - 1/2) 



o 



27ri / r <s ' 7r 

x e O{ W H +W H^ +W H^ +W ) ^ wH + w 3 t 2/3 + w 2 t l/3 + w) 

/l " ^du;|e^ oKw)){1+Xl)+ * 2/ ^ 



o 



7T 

x 0( W 4 i + w¥ /3 + wY /3 +«;) (5.16) 



for some XiiX2,Xz which can be made as small as desired by choosing 
S small enough. At the integration boundary w = 5 the leading term 
is exp(— 8<5 3 t(l + Xi)/3). This easily implies that the integral remains 
bounded as t — > oo. Now we do the change of variable z = t l ^w. The 
t^dw = dz and 0(wH + ... + w) = 0(z A + ... + z)^ 1 ' 3 . The rest of 
the integral is e -3 z3 ( 1 +xi)+nz 2 +c 2 z £ Qr some constants Ci,C2, since the func- 
tion fi do not contains the error terms. The integrand is then t -1 / 3 times 
0(z A + ... + z ) e -h 3 d+xi)+ciz 2 +c 2 z and the integral is on [0,^/3]. The 

e -^z 3 (i+xi) dominates the integral for large z. Thus t 1 ^ 3 x ( I5.16P remains 
finite in the t — > oo limit. Therefore the above estimate of the error term 
becomes F(-l/2)0(t~ 1 / 3 ). 

The final step is to compute z |^ 3 - J T s dvF(v). Extending 5 to oo we only 

make an error of order F{— l/2)(9(e~ M *) for some < fi ~ 5 3 and this leads 
to the integration along the path = {e~ lnssn( - w ^ 3 \w\,w G R}. Therefore 



/ dvF(v) = F(-l/2)—— / dw2e- 8w t/3 

M 2m Jo 



-t 1/3 f ~ _^l/3 e -i7r/3 

2vri 



The change of variable z = 2t 1 ^ 3 e 17T / 3 w leads then to 



27ri r s —Am 

-F(-l/2)0(e^). (5.18) 
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Finally we use an Airy function representation 

[ dve v3/3+av2+bv = Ai(a 2 - b) exp(2a 3 /3 - ab) (5.19) 

~ 2m J loo 

to obtain the final result 

\2\ p 2(« 2 -«i) 3 /3+(si+s 2 )(«2-Mi) 



, lim -WTTT^ I dvF ^ = Ai ^ Sl + S2 + ^ _ Ul ^ e2 
t^co F{-l/2) J To 



(5.20) 
□ 



A Compact form for the extended kernel 

In this Appendix we show that the entries of the compact form of the kernel 
(11.41) agree with (12.101) . Let us introduce some notations. Let Q be the 
multiplication operator by the position, D be the differentiation operator, 
and let A be the Laplacian. On Schwarz test functions / G <S(R), 

(Qf)(x)=xf(x), (Df)(x) = ^f(x), (Af)(x) = ^f(x). (A.l) 

Moreover, denote by K x the operator with kernel K\(x,y) = Ai(x + y + A) 
and the Airy operator Ha = —A + Q. 

We will apply Baker- Campbell-Haussdorf formula. If [A, [A, B]] = ct and 
[B, [A, B]] = c'l for some constant c, c', then 

e A e B = e A ^ B +\[ A M^[A\AM-^[B,{A,B\\ (A.2) 

from which follows, for [A, B] = ct, 

e A+B =e A e B e -±[A,B]_ (A _ 3) 

Moreover, we will use the property of Airy functions 

M"(x + y) = (x + y)Ai(x + y). (A.4) 
We collect some useful properties in the following lemma. 
Lemma A.l. 

1) Commutation relations: [Q, D] = —t, [Q,A] = —2D, 
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2) H A K X = -K X (X1 + Q), 

3) e tD K x = K\ +t and K x e tD = K x „ t , 

4) e~ tA K x = e-l t3 - xt e- t ®K x+t 2e- t Q. 



Proof. 1) By applying to f E «S(R) we get [Q,D]f(x) = —f(x) and 

[Q,A]f(x) = -2f'(x) = -2Df(x). 

2) We apply the definition of Ha and use ( 1A.4I) to get 



(H A K x f)(x) = - J dyAi(x+y+X)(y+\)f(y) = -(K x (\l+Q)f)(x). (A.5) 

3) Follows from (e tD f)(x) = f(x + t). 

4) We use the property 2) and fTA~2l with A = -tA and B = -t(-A + Q), 

to get 

e - tA K x = e' tA e^- A+Q ^e tHA K x = e~ tQ+t2D K x e- tQ e~ M (A.6) 
Then apply f[Ql) with A = -tQ and B = t 2 D to get 4). □ 

What we have to compute explicitly is e~ UlA Koe U2A . From 4) of 
Lemma [A. II we have 
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e~ UlA K e U2A = e-3 u ^e- UlQ K u 2e- UlQ e U2A . (A.7) 
The last part can be rewritten as 

g-Mi<9 e u 2A _ e (u 2 A-u 1 Q-^U2U 2 -u 1 u 2 D)+(2u 1 U2D) 



by ([Oil with A = and 5 = u 2 A. Then using flA,3p with 

A = u 2 A — U\Q — \u 2 u\ — U\U 2 D and B = 1u\U 2 D we obtain 

g-«iQg«2A _ e 2Miu 2 15gM2A e -«ig e «2n? ^ g^ 

Plugging this back into ( 1A.7h we have 

e- UlA K e U2A = e -h" +U2U " e - UlQ K u 2e 2u2UlD e U2A e- UlQ . (A. 10) 

Then we apply 3) of Lemma TA.lt namely K u 2e 2u2UlD = K u 2_ 2u2Ul , and we 
exchange the order of K. and e U2A because are both symmetric and apply 4) 
of Lemma IA.11 This results into 

e- UlA K e U2A = e -l {ui - U2f ''e- (ui - U2)Q K {ui _ U2) 2e- {ui - U2)Q '. (A.ll) 
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Explicitly 

(e- ulA K e U2A )(s u s 2 ) = e l (n2 - Ul)3 e (n2 - ul)(si+S2) Ai( Sl + s 2 + ( Ul - u 2 f). 

(A.12) 

Thus we showed how the second term in HI .41) leads to the corresponding one 
in (11.5j) . It remains the first one, (e^ U2 ~ Ul ^ A )(si, s 2 ), for u 2 > u\. This is just 
the one-dimensional heat kernel, for which it is well known that (see e.g. [26]) 

(e(—)A)( SljS2 ) = / exp • (A.13) 

B Charlier polynomials 

In this Appendix we explain a constructive method to do the orthogonal- 
ization. Let C n (x,t) be the Charlier polynomial of degree n. They are 
orthogonal polynomials with respect to the weight on {0, 1, . . .} given by 

w t (z) = e~H z /z\ (B.l) 

which are traditionally normalized via 

v — v Tl\ 

2_^C n (z,t)C m (z,t)w t (z) = — 5 ntm (B.2) 

V 

z>0 

or, equivalently, C n (z, t) = (—l/t) n z n H — • . They can be expressed in terms 
of hypergeometric functions 

C n {x, t) = 2 F {-n, -x; ; -1/t) (B.3) 

and satisfy the recurrence relation 

~C n {x-l,t) = C n (x,t)-C n+1 (x 1 t). (B.4) 



From the generating function of the Charlier polynomials 

£ ^2^-v n = e v {l - v/t) x (B.5) 



n>0 



one gets the integral representation 



1 „ , , 1 / dve v (l -v/tY ^ , 

C n (z,t) = —<b 1 ' ■ (B.6) 



n\ 2-k'i J Tq v 
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For a good reference on orthogonal polynomials, see [20] . 

It is not too difficult to see that the functions defined in Lemma 14.11 
can be expressed in terms of the Charlier orthogonal polynomials as 

*Z(z) = 0^C k (z-k,t). (B.7) 



Using the recurrence relation ( 1B.4I) repeatedly we obtain 

2k 



^(z)=w t (z)J2Sk,iC l (z,t), (B.8) 

1=0 

where the entries of the matrix S are 

Notice that S is not a square matrix. From this it follows that the polyno- 
mials $^ which satisfy 

'£*S(z)9»(z)=5 kJ , (B.10) 



z>0 



are given by 



N-l 



= (B.n) 

1=0 

where by S^ 1 we mean the (i, j)-entry of the inverse of the square matrix 

S = [Sij]o<i,j<N-i obtained by restricting S to the first N indices. The 
first main difficulty is to obtain the inverse of S. After some work we could 
determine it, namely 

(* /),/.• ™ 

with the identification Sqq = 1 and the convention that the RHS of (IB. 12|) 
is zero when i > j. 

At this point we substitute ( IB. 121) into ( IB. Ill) , perform the summation, 
and finally change the variable v = —wt. The final result is the biorthogonal 
functions $^ reported in (14. 3p . 
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